

%% election outcomes by PRI-PVEM coalition configuration

e12 = readtable('../raw_data/election_results_2012.csv');

% order observations by candidate-coalition choice
d0 = find(e12.candCM==0); N0 = size(d0,1);
d1 = find(e12.candCM==1); N1 = size(d1,1);
d2 = find(e12.candCM==2); N2 = size(d2,1); N = N0+N1+N2;
e12 = e12([d0;d1;d2],:);

votes = [e12.MP,e12.NA,e12.PVEM+0.5*e12.CM,e12.PRI+0.5*e12.CM,e12.PAN];
[~,winner] = max(votes,[],2);
[~,winner(N0+(1:N1))] = max([votes(N0+(1:N1),1:2), ...
    zeros(N1,1),sum(votes(N0+(1:N1),3:4),2),votes(N0+(1:N1),5)],[],2);
[~,winner(N0+N1+(1:N2))] = max([votes(N0+N1+(1:N2),1:2), ...
    sum(votes(N0+N1+(1:N2),3:4),2),zeros(N2,1),votes(N0+N1+(1:N2),5)],[],2);

table2 = table;
table2.party = {'MP';'NA';'PVEM';'PRI';'PAN'};
table2.victory_rate_distinct = round(100 * [sum(winner(1:N0)==1);sum(winner(1:N0)==2); ...
    sum(winner(1:N0)==3);sum(winner(1:N0)==4);sum(winner(1:N0)==5)] / N0,1);
table2.vote_share_distinct = round(100 * mean(votes(1:N0,:) ./ sum(votes(1:N0,:),2),1)',1);
table2.victory_rate_jointPRI = round(100 * [sum(winner(N0+(1:N1))==1);sum(winner(N0+(1:N1))==2); ...
    sum(winner(N0+(1:N1))==3);sum(winner(N0+(1:N1))==4);sum(winner(N0+(1:N1))==5)] / N1,1);
table2.vote_share_jointPRI = round(100 * mean(votes(N0+(1:N1),:) ./ sum(votes(N0+(1:N1),:),2),1)',1);
table2.victory_rate_jointPVEM = round(100 * [sum(winner(N0+N1+(1:N2))==1);sum(winner(N0+N1+(1:N2))==2); ...
    sum(winner(N0+N1+(1:N2))==3);sum(winner(N0+N1+(1:N2))==4);sum(winner(N0+N1+(1:N2))==5)] / N2,1);
table2.vote_share_jointPVEM = round(100 * mean(votes(N0+N1+(1:N2),:) ./ sum(votes(N0+N1+(1:N2),:),2),1)',1);

disp(' ')
disp('for Table 2:')
disp(table2)


%% votes in support of PRI-PVEM coalition candidates

table3 = table;
table3.PRvote = {'PVEM';'PRI';'50-50 Split'};
table3.jointPRI = round(100 * mean([e12.PVEM(N0+(1:N1)),e12.PRI(N0+(1:N1)),e12.CM(N0+(1:N1))] ./ ...
    sum(votes(N0+(1:N1),:),2),1)',1);
table3.jointPVEM = round(100 * mean([e12.PVEM(N0+N1+(1:N2)),e12.PRI(N0+N1+(1:N2)),e12.CM(N0+N1+(1:N2))] ./ ...
    sum(votes(N0+N1+(1:N2),:),2),1)',1);

disp(' ')
disp('for Table 3:')
disp(table3)


%% distribution of PRI-PVEM vote split

% 50-50 split
histogram(100 * e12.CM(N0+(1:N1)) ./ ...
    (e12.PRI(N0+(1:N1)) + e12.PVEM(N0+(1:N1)) + e12.CM(N0+(1:N1))), ...
    7,'Normalization','probability')
hold on
histogram(100 * e12.CM(N0+N1+(1:N2)) ./ ...
    (e12.PRI(N0+N1+(1:N2)) + e12.PVEM(N0+N1+(1:N2)) + e12.CM(N0+N1+(1:N2))), ...
    7,'Normalization','probability')
ax = gca;
title('50-50 PR Split Among PRI-PVEM Coalition Supporters','FontSize',16,'Interpreter','latex')
legend({'PRI Joint Candidate','PVEM Joint Candidate'},'FontSize',14,'Interpreter','latex')
ax.YLabel.String = 'Relative Frequency';
ax.YLabel.FontSize = 14;
ax.YLabel.Interpreter = 'latex';
ax.YLim = [0,0.35];
ax.XLabel.String = 'Percentage of Coalition Candidate''s Votes';
ax.XLabel.FontSize = 14;
ax.XLabel.Interpreter = 'latex';
saveas(gcf,'../figures/figA2c.png')
hold off

% PRI vote
histogram(100 * e12.PRI(N0+(1:N1)) ./ ...
    (e12.PRI(N0+(1:N1)) + e12.PVEM(N0+(1:N1)) + e12.CM(N0+(1:N1))), ...
    8,'Normalization','probability')
hold on
histogram(100 * e12.PRI(N0+N1+(1:N2)) ./ ...
    (e12.PRI(N0+N1+(1:N2)) + e12.PVEM(N0+N1+(1:N2)) + e12.CM(N0+N1+(1:N2))), ...
    8,'Normalization','probability')
ax = gca;
title('PR Vote for PRI Among PRI-PVEM Coalition Supporters','FontSize',16,'Interpreter','latex')
legend({'PRI Joint Candidate','PVEM Joint Candidate'},'FontSize',14,'Interpreter','latex')
ax.YLabel.String = 'Relative Frequency';
ax.YLabel.FontSize = 14;
ax.YLabel.Interpreter = 'latex';
ax.YLim = [0,0.475];
ax.XLabel.String = 'Percentage of Coalition Candidate''s Votes';
ax.XLabel.FontSize = 14;
ax.XLabel.Interpreter = 'latex';
saveas(gcf,'../figures/figA2a.png')
hold off

% PVEM vote
histogram(100 * e12.PVEM(N0+(1:N1)) ./ ...
    (e12.PRI(N0+(1:N1)) + e12.PVEM(N0+(1:N1)) + e12.CM(N0+(1:N1))), ...
    9,'Normalization','probability')
hold on
histogram(100 * e12.PVEM(N0+N1+(1:N2)) ./ ...
    (e12.PRI(N0+N1+(1:N2)) + e12.PVEM(N0+N1+(1:N2)) + e12.CM(N0+N1+(1:N2))), ...
    9,'Normalization','probability')
ax = gca;
title('PR Vote for PVEM Among PRI-PVEM Coalition Supporters','FontSize',16,'Interpreter','latex')
legend({'PRI Joint Candidate','PVEM Joint Candidate'},'FontSize',14,'Interpreter','latex')
ax.YLabel.String = 'Relative Frequency';
ax.YLabel.FontSize = 14;
ax.YLabel.Interpreter = 'latex';
ax.XLabel.String = 'Percentage of Coalition Candidate''s Votes';
ax.XLabel.FontSize = 14;
ax.XLabel.Interpreter = 'latex';
saveas(gcf,'../figures/figA2b.png')
hold off


%% district characteristics

dChar = readtable('../raw_data/district_characteristics.csv');
dChar = dChar([d0;d1;d2],:);

tableA1 = table;
tableA1.demographics = {'Female Head of Household (%)';'Pop. over 60 (% Voting-Age Pop.)';'Rural Neighborhoods (%)'};
tableA1.mean_distinct = round(100 * [mean(dChar.femaleHH(1:N0)); ...
    mean(dChar.over60(1:N0));mean(dChar.rural(1:N0))],1);
tableA1.std_distinct = round(100 * [std(dChar.femaleHH(1:N0)); ...
    std(dChar.over60(1:N0));std(dChar.rural(1:N0))],1);
tableA1.mean_jointPRI = round(100 * [mean(dChar.femaleHH(N0+(1:N1))); ...
    mean(dChar.over60(N0+(1:N1)));mean(dChar.rural(N0+(1:N1)))],1);
tableA1.std_jointPRI = round(100 * [std(dChar.femaleHH(N0+(1:N1))); ...
    std(dChar.over60(N0+(1:N1)));std(dChar.rural(N0+(1:N1)))],1);
tableA1.mean_jointPVEM = round(100 * [mean(dChar.femaleHH(N0+N1+(1:N2))); ...
    mean(dChar.over60(N0+N1+(1:N2)));mean(dChar.rural(N0+N1+(1:N2)))],1);
tableA1.std_jointPVEM = round(100 * [std(dChar.femaleHH(N0+N1+(1:N2))); ...
    std(dChar.over60(N0+N1+(1:N2)));std(dChar.rural(N0+N1+(1:N2)))],1);

disp(' ')
disp('for Table A1:')
disp(tableA1)


%% 2012 campaing expenditures

cs12 = readtable('../raw_data/campaign_spending_2012.csv');
cs12 = cs12([d0;d1;d2],:);

tableA2 = table;
tableA2.party = {'MP';'NA';'PVEM';'PRI';'PAN'};
tableA2.mean_distinct = round(10 * [mean(cs12.MP(1:N0));mean(cs12.NA(1:N0)); ...
    mean(cs12.PVEM(1:N0));mean(cs12.PRI(1:N0));mean(cs12.PAN(1:N0))],1);
tableA2.std_distinct = round(10 * [std(cs12.MP(1:N0));std(cs12.NA(1:N0)); ...
    std(cs12.PVEM(1:N0));std(cs12.PRI(1:N0));std(cs12.PAN(1:N0))],1);
tableA2.mean_jointPRI = round(10 * [mean(cs12.MP(N0+(1:N1)));mean(cs12.NA(N0+(1:N1))); ...
    mean(cs12.PVEM(N0+(1:N1)));mean(cs12.CM(N0+(1:N1)));mean(cs12.PAN(N0+(1:N1)))],1);
tableA2.std_jointPRI = round(10 * [std(cs12.MP(N0+(1:N1)));std(cs12.NA(N0+(1:N1))); ...
    std(cs12.PVEM(N0+(1:N1)));std(cs12.CM(N0+(1:N1)));std(cs12.PAN(N0+(1:N1)))],1);
tableA2.mean_jointPVEM = round(10 * [mean(cs12.MP(N0+N1+(1:N2)));mean(cs12.NA(N0+N1+(1:N2))); ...
    mean(cs12.CM(N0+N1+(1:N2)));mean(cs12.PRI(N0+N1+(1:N2)));mean(cs12.PAN(N0+N1+(1:N2)))],1);
tableA2.std_jointPVEM = round(10 * [std(cs12.MP(N0+N1+(1:N2)));std(cs12.NA(N0+N1+(1:N2))); ...
    std(cs12.CM(N0+N1+(1:N2)));std(cs12.PRI(N0+N1+(1:N2)));std(cs12.PAN(N0+N1+(1:N2)))],1);

disp(' ')
disp('for Table A2:')
disp(tableA2)


%% previous candidate experience

warning('off','all') % suppress Matlab table name format warnings
FPTP2003 = readtable('../raw_data/2003FPTP.xlsx');
FPTP2003 = FPTP2003(strcmp(FPTP2003.CALIDAD_CANDIDATO,'PROPIETARIO'),:);
PR2003 = readtable('../raw_data/2003PR.xlsx');
PR2003.Properties.VariableNames(4:end) = FPTP2003.Properties.VariableNames(4:end);
PR2003 = PR2003(strcmp(PR2003.CALIDAD_CANDIDATO,'PROPIETARIO'),:);
FPTP2006 = readtable('../raw_data/2006FPTP.xlsx');
FPTP2006.Properties.VariableNames(4:end) = FPTP2003.Properties.VariableNames(4:end);
FPTP2006 = FPTP2006(strcmp(FPTP2006.CALIDAD_CANDIDATO,'PROPIETARIO'),:);
PR2006 = readtable('../raw_data/2006PR.xlsx');
PR2006.Properties.VariableNames(4:end) = FPTP2003.Properties.VariableNames(4:end);
PR2006 = PR2006(strcmp(PR2006.CALIDAD_CANDIDATO,'PROPIETARIO'),:);
FPTP2006S = readtable('../raw_data/2006FPTPsenate.xlsx');
FPTP2006S.Properties.VariableNames(4:end) = FPTP2003.Properties.VariableNames(4:end);
FPTP2006S = FPTP2006S(strcmp(FPTP2006S.CALIDAD_CANDIDATO,'PROPIETARIO'),:);
PR2006S = readtable('../raw_data/2006PRsenate.xlsx');
PR2006S.Properties.VariableNames(3:end) = FPTP2003.Properties.VariableNames(4:end);
PR2006S = PR2006S(strcmp(PR2006S.CALIDAD_CANDIDATO,'PROPIETARIO'),:);
FPTP2009 = readtable('../raw_data/2009FPTP.xlsx');
FPTP2009 = FPTP2009(:,[1:4,6,7,5]);
FPTP2009.Properties.VariableNames(4:end) = FPTP2003.Properties.VariableNames(4:end);
FPTP2009 = FPTP2009(strcmp(FPTP2009.CALIDAD_CANDIDATO,'PROPIETARIO'),:);
PR2009 = readtable('../raw_data/2009PR.xlsx');
PR2009 = PR2009(:,[1,2,4,3,6,7,5]);
PR2009.Properties.VariableNames(4:end) = FPTP2003.Properties.VariableNames(4:end);
PR2009 = PR2009(strcmp(PR2009.CALIDAD_CANDIDATO,'PROPIETARIO'),:);
FPTP2012 = readtable('../raw_data/2012FPTP.xlsx');
FPTP2012 = FPTP2012(:,[1:4,6,7,5]);
FPTP2012.Properties.VariableNames(4:end) = FPTP2003.Properties.VariableNames(4:end);
FPTP2012 = FPTP2012(strcmp(FPTP2012.CALIDAD_CANDIDATO,'PROPIETARIO'),:);
PR2012 = readtable('../raw_data/2012PR.xlsx');
PR2012 = PR2012(:,[1:4,6,7,5]);
PR2012.Properties.VariableNames(4:end) = FPTP2003.Properties.VariableNames(4:end);
PR2012 = PR2012(strcmp(PR2012.CALIDAD_CANDIDATO,'PROPIETARIO'),:);

FPTP2012.ran2009 = zeros(size(FPTP2012,1),1);
FPTP2012.ran2006 = FPTP2012.ran2009;
FPTP2012.ran2003 = FPTP2012.ran2009;
PR2012.ran2009 = zeros(size(PR2012,1),1);
PR2012.ran2006 = PR2012.ran2009;
PR2012.ran2003 = PR2012.ran2009;
MP_FPTP = FPTP2012(strcmp(FPTP2012.PARTIDOOCOALICI_N,'MOVIMIENTO PROGRESISTA'),:);
MP_PR = PR2012(strcmp(PR2012.PARTIDO,'PRD')|strcmp(PR2012.PARTIDO,'PT')|strcmp(PR2012.PARTIDO,'MOVIMIENTO CIUDADANO'),:);
NA_FPTP = FPTP2012(strcmp(FPTP2012.PARTIDOOCOALICI_N,'NUEVA ALIANZA'),:);
NA_PR = PR2012(strcmp(PR2012.PARTIDO,'NUEVA ALIANZA'),:);
PVEM_FPTP = FPTP2012(strcmp(FPTP2012.PARTIDOOCOALICI_N,'PVEM'),:);
PVEM_PR = PR2012(strcmp(PR2012.PARTIDO,'PVEM'),:);
PRI_FPTP = FPTP2012(strcmp(FPTP2012.PARTIDOOCOALICI_N,'PRI'),:);
PRI_PR = PR2012(strcmp(PR2012.PARTIDO,'PRI'),:);
CM_FPTP = FPTP2012(strcmp(FPTP2012.PARTIDOOCOALICI_N,'COMPROMISO POR MÉXICO'),:);
CM_affil = readtable('../raw_data/election_results_2012.csv');
CM_affil = CM_affil.candCM(CM_affil.candCM>0);
PVEM_FPTP = [PVEM_FPTP;CM_FPTP(CM_affil==2,:)];
PRI_FPTP = [PRI_FPTP;CM_FPTP(CM_affil==1,:)];
clear CM_FPTP CM_affil
PAN_FPTP = FPTP2012(strcmp(FPTP2012.PARTIDOOCOALICI_N,'PAN'),:);
PAN_PR = PR2012(strcmp(PR2012.PARTIDO,'PAN'),:);

% MP
for i = 1:size(MP_FPTP,1)
    MP_FPTP.ran2009(i) = MP_FPTP.ran2009(i)+...
        sum(strcmp(MP_FPTP.APELLIDO_PATERNO{i},FPTP2009.APELLIDO_PATERNO)&...
        strcmp(MP_FPTP.APELLIDO_MATERNO{i},FPTP2009.APELLIDO_MATERNO)&...
        strcmp(MP_FPTP.NOMBRE{i},FPTP2009.NOMBRE))+...
        sum(strcmp(MP_FPTP.APELLIDO_PATERNO{i},PR2009.APELLIDO_PATERNO)&...
        strcmp(MP_FPTP.APELLIDO_MATERNO{i},PR2009.APELLIDO_MATERNO)&...
        strcmp(MP_FPTP.NOMBRE{i},PR2009.NOMBRE));
    MP_FPTP.ran2006(i) = MP_FPTP.ran2006(i)+...
        sum(strcmp(MP_FPTP.APELLIDO_PATERNO{i},FPTP2006.APELLIDO_PATERNO)&...
        strcmp(MP_FPTP.APELLIDO_MATERNO{i},FPTP2006.APELLIDO_MATERNO)&...
        strcmp(MP_FPTP.NOMBRE{i},FPTP2006.NOMBRE))+...
        sum(strcmp(MP_FPTP.APELLIDO_PATERNO{i},PR2006.APELLIDO_PATERNO)&...
        strcmp(MP_FPTP.APELLIDO_MATERNO{i},PR2006.APELLIDO_MATERNO)&...
        strcmp(MP_FPTP.NOMBRE{i},PR2006.NOMBRE))+...
        sum(strcmp(MP_FPTP.APELLIDO_PATERNO{i},FPTP2006S.APELLIDO_PATERNO)&...
        strcmp(MP_FPTP.APELLIDO_MATERNO{i},FPTP2006S.APELLIDO_MATERNO)&...
        strcmp(MP_FPTP.NOMBRE{i},FPTP2006S.NOMBRE))+...
        sum(strcmp(MP_FPTP.APELLIDO_PATERNO{i},PR2006S.APELLIDO_PATERNO)&...
        strcmp(MP_FPTP.APELLIDO_MATERNO{i},PR2006S.APELLIDO_MATERNO)&...
        strcmp(MP_FPTP.NOMBRE{i},PR2006S.NOMBRE));
    MP_FPTP.ran2003(i) = MP_FPTP.ran2003(i)+...
        sum(strcmp(MP_FPTP.APELLIDO_PATERNO{i},FPTP2003.APELLIDO_PATERNO)&...
        strcmp(MP_FPTP.APELLIDO_MATERNO{i},FPTP2003.APELLIDO_MATERNO)&...
        strcmp(MP_FPTP.NOMBRE{i},FPTP2003.NOMBRE))+...
        sum(strcmp(MP_FPTP.APELLIDO_PATERNO{i},PR2003.APELLIDO_PATERNO)&...
        strcmp(MP_FPTP.APELLIDO_MATERNO{i},PR2003.APELLIDO_MATERNO)&...
        strcmp(MP_FPTP.NOMBRE{i},PR2003.NOMBRE));
end
for i = 1:size(MP_PR,1)
    MP_PR.ran2009(i) = MP_PR.ran2009(i)+...
        sum(strcmp(MP_PR.APELLIDO_PATERNO{i},FPTP2009.APELLIDO_PATERNO)&...
        strcmp(MP_PR.APELLIDO_MATERNO{i},FPTP2009.APELLIDO_MATERNO)&...
        strcmp(MP_PR.NOMBRE{i},FPTP2009.NOMBRE))+...
        sum(strcmp(MP_PR.APELLIDO_PATERNO{i},PR2009.APELLIDO_PATERNO)&...
        strcmp(MP_PR.APELLIDO_MATERNO{i},PR2009.APELLIDO_MATERNO)&...
        strcmp(MP_PR.NOMBRE{i},PR2009.NOMBRE));
    MP_PR.ran2006(i) = MP_PR.ran2006(i)+...
        sum(strcmp(MP_PR.APELLIDO_PATERNO{i},FPTP2006.APELLIDO_PATERNO)&...
        strcmp(MP_PR.APELLIDO_MATERNO{i},FPTP2006.APELLIDO_MATERNO)&...
        strcmp(MP_PR.NOMBRE{i},FPTP2006.NOMBRE))+...
        sum(strcmp(MP_PR.APELLIDO_PATERNO{i},PR2006.APELLIDO_PATERNO)&...
        strcmp(MP_PR.APELLIDO_MATERNO{i},PR2006.APELLIDO_MATERNO)&...
        strcmp(MP_PR.NOMBRE{i},PR2006.NOMBRE))+...
        sum(strcmp(MP_PR.APELLIDO_PATERNO{i},FPTP2006S.APELLIDO_PATERNO)&...
        strcmp(MP_PR.APELLIDO_MATERNO{i},FPTP2006S.APELLIDO_MATERNO)&...
        strcmp(MP_PR.NOMBRE{i},FPTP2006S.NOMBRE))+...
        sum(strcmp(MP_PR.APELLIDO_PATERNO{i},PR2006S.APELLIDO_PATERNO)&...
        strcmp(MP_PR.APELLIDO_MATERNO{i},PR2006S.APELLIDO_MATERNO)&...
        strcmp(MP_PR.NOMBRE{i},PR2006S.NOMBRE));
    MP_PR.ran2003(i) = MP_PR.ran2003(i)+...
        sum(strcmp(MP_PR.APELLIDO_PATERNO{i},FPTP2003.APELLIDO_PATERNO)&...
        strcmp(MP_PR.APELLIDO_MATERNO{i},FPTP2003.APELLIDO_MATERNO)&...
        strcmp(MP_PR.NOMBRE{i},FPTP2003.NOMBRE))+...
        sum(strcmp(MP_PR.APELLIDO_PATERNO{i},PR2003.APELLIDO_PATERNO)&...
        strcmp(MP_PR.APELLIDO_MATERNO{i},PR2003.APELLIDO_MATERNO)&...
        strcmp(MP_PR.NOMBRE{i},PR2003.NOMBRE));
end

% NA
for i = 1:size(NA_FPTP,1)
    NA_FPTP.ran2009(i) = NA_FPTP.ran2009(i)+...
        sum(strcmp(NA_FPTP.APELLIDO_PATERNO{i},FPTP2009.APELLIDO_PATERNO)&...
        strcmp(NA_FPTP.APELLIDO_MATERNO{i},FPTP2009.APELLIDO_MATERNO)&...
        strcmp(NA_FPTP.NOMBRE{i},FPTP2009.NOMBRE))+...
        sum(strcmp(NA_FPTP.APELLIDO_PATERNO{i},PR2009.APELLIDO_PATERNO)&...
        strcmp(NA_FPTP.APELLIDO_MATERNO{i},PR2009.APELLIDO_MATERNO)&...
        strcmp(NA_FPTP.NOMBRE{i},PR2009.NOMBRE));
    NA_FPTP.ran2006(i) = NA_FPTP.ran2006(i)+...
        sum(strcmp(NA_FPTP.APELLIDO_PATERNO{i},FPTP2006.APELLIDO_PATERNO)&...
        strcmp(NA_FPTP.APELLIDO_MATERNO{i},FPTP2006.APELLIDO_MATERNO)&...
        strcmp(NA_FPTP.NOMBRE{i},FPTP2006.NOMBRE))+...
        sum(strcmp(NA_FPTP.APELLIDO_PATERNO{i},PR2006.APELLIDO_PATERNO)&...
        strcmp(NA_FPTP.APELLIDO_MATERNO{i},PR2006.APELLIDO_MATERNO)&...
        strcmp(NA_FPTP.NOMBRE{i},PR2006.NOMBRE))+...
        sum(strcmp(NA_FPTP.APELLIDO_PATERNO{i},FPTP2006S.APELLIDO_PATERNO)&...
        strcmp(NA_FPTP.APELLIDO_MATERNO{i},FPTP2006S.APELLIDO_MATERNO)&...
        strcmp(NA_FPTP.NOMBRE{i},FPTP2006S.NOMBRE))+...
        sum(strcmp(NA_FPTP.APELLIDO_PATERNO{i},PR2006S.APELLIDO_PATERNO)&...
        strcmp(NA_FPTP.APELLIDO_MATERNO{i},PR2006S.APELLIDO_MATERNO)&...
        strcmp(NA_FPTP.NOMBRE{i},PR2006S.NOMBRE));
    NA_FPTP.ran2003(i) = NA_FPTP.ran2003(i)+...
        sum(strcmp(NA_FPTP.APELLIDO_PATERNO{i},FPTP2003.APELLIDO_PATERNO)&...
        strcmp(NA_FPTP.APELLIDO_MATERNO{i},FPTP2003.APELLIDO_MATERNO)&...
        strcmp(NA_FPTP.NOMBRE{i},FPTP2003.NOMBRE))+...
        sum(strcmp(NA_FPTP.APELLIDO_PATERNO{i},PR2003.APELLIDO_PATERNO)&...
        strcmp(NA_FPTP.APELLIDO_MATERNO{i},PR2003.APELLIDO_MATERNO)&...
        strcmp(NA_FPTP.NOMBRE{i},PR2003.NOMBRE));
end
for i = 1:size(NA_PR,1)
    NA_PR.ran2009(i) = NA_PR.ran2009(i)+...
        sum(strcmp(NA_PR.APELLIDO_PATERNO{i},FPTP2009.APELLIDO_PATERNO)&...
        strcmp(NA_PR.APELLIDO_MATERNO{i},FPTP2009.APELLIDO_MATERNO)&...
        strcmp(NA_PR.NOMBRE{i},FPTP2009.NOMBRE))+...
        sum(strcmp(NA_PR.APELLIDO_PATERNO{i},PR2009.APELLIDO_PATERNO)&...
        strcmp(NA_PR.APELLIDO_MATERNO{i},PR2009.APELLIDO_MATERNO)&...
        strcmp(NA_PR.NOMBRE{i},PR2009.NOMBRE));
    NA_PR.ran2006(i) = NA_PR.ran2006(i)+...
        sum(strcmp(NA_PR.APELLIDO_PATERNO{i},FPTP2006.APELLIDO_PATERNO)&...
        strcmp(NA_PR.APELLIDO_MATERNO{i},FPTP2006.APELLIDO_MATERNO)&...
        strcmp(NA_PR.NOMBRE{i},FPTP2006.NOMBRE))+...
        sum(strcmp(NA_PR.APELLIDO_PATERNO{i},PR2006.APELLIDO_PATERNO)&...
        strcmp(NA_PR.APELLIDO_MATERNO{i},PR2006.APELLIDO_MATERNO)&...
        strcmp(NA_PR.NOMBRE{i},PR2006.NOMBRE))+...
        sum(strcmp(NA_PR.APELLIDO_PATERNO{i},FPTP2006S.APELLIDO_PATERNO)&...
        strcmp(NA_PR.APELLIDO_MATERNO{i},FPTP2006S.APELLIDO_MATERNO)&...
        strcmp(NA_PR.NOMBRE{i},FPTP2006S.NOMBRE))+...
        sum(strcmp(NA_PR.APELLIDO_PATERNO{i},PR2006S.APELLIDO_PATERNO)&...
        strcmp(NA_PR.APELLIDO_MATERNO{i},PR2006S.APELLIDO_MATERNO)&...
        strcmp(NA_PR.NOMBRE{i},PR2006S.NOMBRE));
    NA_PR.ran2003(i) = NA_PR.ran2003(i)+...
        sum(strcmp(NA_PR.APELLIDO_PATERNO{i},FPTP2003.APELLIDO_PATERNO)&...
        strcmp(NA_PR.APELLIDO_MATERNO{i},FPTP2003.APELLIDO_MATERNO)&...
        strcmp(NA_PR.NOMBRE{i},FPTP2003.NOMBRE))+...
        sum(strcmp(NA_PR.APELLIDO_PATERNO{i},PR2003.APELLIDO_PATERNO)&...
        strcmp(NA_PR.APELLIDO_MATERNO{i},PR2003.APELLIDO_MATERNO)&...
        strcmp(NA_PR.NOMBRE{i},PR2003.NOMBRE));
end

% PVEM
for i = 1:size(PVEM_FPTP,1)
    PVEM_FPTP.ran2009(i) = PVEM_FPTP.ran2009(i)+...
        sum(strcmp(PVEM_FPTP.APELLIDO_PATERNO{i},FPTP2009.APELLIDO_PATERNO)&...
        strcmp(PVEM_FPTP.APELLIDO_MATERNO{i},FPTP2009.APELLIDO_MATERNO)&...
        strcmp(PVEM_FPTP.NOMBRE{i},FPTP2009.NOMBRE))+...
        sum(strcmp(PVEM_FPTP.APELLIDO_PATERNO{i},PR2009.APELLIDO_PATERNO)&...
        strcmp(PVEM_FPTP.APELLIDO_MATERNO{i},PR2009.APELLIDO_MATERNO)&...
        strcmp(PVEM_FPTP.NOMBRE{i},PR2009.NOMBRE));
    PVEM_FPTP.ran2006(i) = PVEM_FPTP.ran2006(i)+...
        sum(strcmp(PVEM_FPTP.APELLIDO_PATERNO{i},FPTP2006.APELLIDO_PATERNO)&...
        strcmp(PVEM_FPTP.APELLIDO_MATERNO{i},FPTP2006.APELLIDO_MATERNO)&...
        strcmp(PVEM_FPTP.NOMBRE{i},FPTP2006.NOMBRE))+...
        sum(strcmp(PVEM_FPTP.APELLIDO_PATERNO{i},PR2006.APELLIDO_PATERNO)&...
        strcmp(PVEM_FPTP.APELLIDO_MATERNO{i},PR2006.APELLIDO_MATERNO)&...
        strcmp(PVEM_FPTP.NOMBRE{i},PR2006.NOMBRE))+...
        sum(strcmp(PVEM_FPTP.APELLIDO_PATERNO{i},FPTP2006S.APELLIDO_PATERNO)&...
        strcmp(PVEM_FPTP.APELLIDO_MATERNO{i},FPTP2006S.APELLIDO_MATERNO)&...
        strcmp(PVEM_FPTP.NOMBRE{i},FPTP2006S.NOMBRE))+...
        sum(strcmp(PVEM_FPTP.APELLIDO_PATERNO{i},PR2006S.APELLIDO_PATERNO)&...
        strcmp(PVEM_FPTP.APELLIDO_MATERNO{i},PR2006S.APELLIDO_MATERNO)&...
        strcmp(PVEM_FPTP.NOMBRE{i},PR2006S.NOMBRE));
    PVEM_FPTP.ran2003(i) = PVEM_FPTP.ran2003(i)+...
        sum(strcmp(PVEM_FPTP.APELLIDO_PATERNO{i},FPTP2003.APELLIDO_PATERNO)&...
        strcmp(PVEM_FPTP.APELLIDO_MATERNO{i},FPTP2003.APELLIDO_MATERNO)&...
        strcmp(PVEM_FPTP.NOMBRE{i},FPTP2003.NOMBRE))+...
        sum(strcmp(PVEM_FPTP.APELLIDO_PATERNO{i},PR2003.APELLIDO_PATERNO)&...
        strcmp(PVEM_FPTP.APELLIDO_MATERNO{i},PR2003.APELLIDO_MATERNO)&...
        strcmp(PVEM_FPTP.NOMBRE{i},PR2003.NOMBRE));
end
for i = 1:size(PVEM_PR,1)
    PVEM_PR.ran2009(i) = PVEM_PR.ran2009(i)+...
        sum(strcmp(PVEM_PR.APELLIDO_PATERNO{i},FPTP2009.APELLIDO_PATERNO)&...
        strcmp(PVEM_PR.APELLIDO_MATERNO{i},FPTP2009.APELLIDO_MATERNO)&...
        strcmp(PVEM_PR.NOMBRE{i},FPTP2009.NOMBRE))+...
        sum(strcmp(PVEM_PR.APELLIDO_PATERNO{i},PR2009.APELLIDO_PATERNO)&...
        strcmp(PVEM_PR.APELLIDO_MATERNO{i},PR2009.APELLIDO_MATERNO)&...
        strcmp(PVEM_PR.NOMBRE{i},PR2009.NOMBRE));
    PVEM_PR.ran2006(i) = PVEM_PR.ran2006(i)+...
        sum(strcmp(PVEM_PR.APELLIDO_PATERNO{i},FPTP2006.APELLIDO_PATERNO)&...
        strcmp(PVEM_PR.APELLIDO_MATERNO{i},FPTP2006.APELLIDO_MATERNO)&...
        strcmp(PVEM_PR.NOMBRE{i},FPTP2006.NOMBRE))+...
        sum(strcmp(PVEM_PR.APELLIDO_PATERNO{i},PR2006.APELLIDO_PATERNO)&...
        strcmp(PVEM_PR.APELLIDO_MATERNO{i},PR2006.APELLIDO_MATERNO)&...
        strcmp(PVEM_PR.NOMBRE{i},PR2006.NOMBRE))+...
        sum(strcmp(PVEM_PR.APELLIDO_PATERNO{i},FPTP2006S.APELLIDO_PATERNO)&...
        strcmp(PVEM_PR.APELLIDO_MATERNO{i},FPTP2006S.APELLIDO_MATERNO)&...
        strcmp(PVEM_PR.NOMBRE{i},FPTP2006S.NOMBRE))+...
        sum(strcmp(PVEM_PR.APELLIDO_PATERNO{i},PR2006S.APELLIDO_PATERNO)&...
        strcmp(PVEM_PR.APELLIDO_MATERNO{i},PR2006S.APELLIDO_MATERNO)&...
        strcmp(PVEM_PR.NOMBRE{i},PR2006S.NOMBRE));
    PVEM_PR.ran2003(i) = PVEM_PR.ran2003(i)+...
        sum(strcmp(PVEM_PR.APELLIDO_PATERNO{i},FPTP2003.APELLIDO_PATERNO)&...
        strcmp(PVEM_PR.APELLIDO_MATERNO{i},FPTP2003.APELLIDO_MATERNO)&...
        strcmp(PVEM_PR.NOMBRE{i},FPTP2003.NOMBRE))+...
        sum(strcmp(PVEM_PR.APELLIDO_PATERNO{i},PR2003.APELLIDO_PATERNO)&...
        strcmp(PVEM_PR.APELLIDO_MATERNO{i},PR2003.APELLIDO_MATERNO)&...
        strcmp(PVEM_PR.NOMBRE{i},PR2003.NOMBRE));
end

% PRI
for i = 1:size(PRI_FPTP,1)
    PRI_FPTP.ran2009(i) = PRI_FPTP.ran2009(i)+...
        sum(strcmp(PRI_FPTP.APELLIDO_PATERNO{i},FPTP2009.APELLIDO_PATERNO)&...
        strcmp(PRI_FPTP.APELLIDO_MATERNO{i},FPTP2009.APELLIDO_MATERNO)&...
        strcmp(PRI_FPTP.NOMBRE{i},FPTP2009.NOMBRE))+...
        sum(strcmp(PRI_FPTP.APELLIDO_PATERNO{i},PR2009.APELLIDO_PATERNO)&...
        strcmp(PRI_FPTP.APELLIDO_MATERNO{i},PR2009.APELLIDO_MATERNO)&...
        strcmp(PRI_FPTP.NOMBRE{i},PR2009.NOMBRE));
    PRI_FPTP.ran2006(i) = PRI_FPTP.ran2006(i)+...
        sum(strcmp(PRI_FPTP.APELLIDO_PATERNO{i},FPTP2006.APELLIDO_PATERNO)&...
        strcmp(PRI_FPTP.APELLIDO_MATERNO{i},FPTP2006.APELLIDO_MATERNO)&...
        strcmp(PRI_FPTP.NOMBRE{i},FPTP2006.NOMBRE))+...
        sum(strcmp(PRI_FPTP.APELLIDO_PATERNO{i},PR2006.APELLIDO_PATERNO)&...
        strcmp(PRI_FPTP.APELLIDO_MATERNO{i},PR2006.APELLIDO_MATERNO)&...
        strcmp(PRI_FPTP.NOMBRE{i},PR2006.NOMBRE))+...
        sum(strcmp(PRI_FPTP.APELLIDO_PATERNO{i},FPTP2006S.APELLIDO_PATERNO)&...
        strcmp(PRI_FPTP.APELLIDO_MATERNO{i},FPTP2006S.APELLIDO_MATERNO)&...
        strcmp(PRI_FPTP.NOMBRE{i},FPTP2006S.NOMBRE))+...
        sum(strcmp(PRI_FPTP.APELLIDO_PATERNO{i},PR2006S.APELLIDO_PATERNO)&...
        strcmp(PRI_FPTP.APELLIDO_MATERNO{i},PR2006S.APELLIDO_MATERNO)&...
        strcmp(PRI_FPTP.NOMBRE{i},PR2006S.NOMBRE));
    PRI_FPTP.ran2003(i) = PRI_FPTP.ran2003(i)+...
        sum(strcmp(PRI_FPTP.APELLIDO_PATERNO{i},FPTP2003.APELLIDO_PATERNO)&...
        strcmp(PRI_FPTP.APELLIDO_MATERNO{i},FPTP2003.APELLIDO_MATERNO)&...
        strcmp(PRI_FPTP.NOMBRE{i},FPTP2003.NOMBRE))+...
        sum(strcmp(PRI_FPTP.APELLIDO_PATERNO{i},PR2003.APELLIDO_PATERNO)&...
        strcmp(PRI_FPTP.APELLIDO_MATERNO{i},PR2003.APELLIDO_MATERNO)&...
        strcmp(PRI_FPTP.NOMBRE{i},PR2003.NOMBRE));
end
for i = 1:size(PRI_PR,1)
    PRI_PR.ran2009(i) = PRI_PR.ran2009(i)+...
        sum(strcmp(PRI_PR.APELLIDO_PATERNO{i},FPTP2009.APELLIDO_PATERNO)&...
        strcmp(PRI_PR.APELLIDO_MATERNO{i},FPTP2009.APELLIDO_MATERNO)&...
        strcmp(PRI_PR.NOMBRE{i},FPTP2009.NOMBRE))+...
        sum(strcmp(PRI_PR.APELLIDO_PATERNO{i},PR2009.APELLIDO_PATERNO)&...
        strcmp(PRI_PR.APELLIDO_MATERNO{i},PR2009.APELLIDO_MATERNO)&...
        strcmp(PRI_PR.NOMBRE{i},PR2009.NOMBRE));
    PRI_PR.ran2006(i) = PRI_PR.ran2006(i)+...
        sum(strcmp(PRI_PR.APELLIDO_PATERNO{i},FPTP2006.APELLIDO_PATERNO)&...
        strcmp(PRI_PR.APELLIDO_MATERNO{i},FPTP2006.APELLIDO_MATERNO)&...
        strcmp(PRI_PR.NOMBRE{i},FPTP2006.NOMBRE))+...
        sum(strcmp(PRI_PR.APELLIDO_PATERNO{i},PR2006.APELLIDO_PATERNO)&...
        strcmp(PRI_PR.APELLIDO_MATERNO{i},PR2006.APELLIDO_MATERNO)&...
        strcmp(PRI_PR.NOMBRE{i},PR2006.NOMBRE))+...
        sum(strcmp(PRI_PR.APELLIDO_PATERNO{i},FPTP2006S.APELLIDO_PATERNO)&...
        strcmp(PRI_PR.APELLIDO_MATERNO{i},FPTP2006S.APELLIDO_MATERNO)&...
        strcmp(PRI_PR.NOMBRE{i},FPTP2006S.NOMBRE))+...
        sum(strcmp(PRI_PR.APELLIDO_PATERNO{i},PR2006S.APELLIDO_PATERNO)&...
        strcmp(PRI_PR.APELLIDO_MATERNO{i},PR2006S.APELLIDO_MATERNO)&...
        strcmp(PRI_PR.NOMBRE{i},PR2006S.NOMBRE));
    PRI_PR.ran2003(i) = PRI_PR.ran2003(i)+...
        sum(strcmp(PRI_PR.APELLIDO_PATERNO{i},FPTP2003.APELLIDO_PATERNO)&...
        strcmp(PRI_PR.APELLIDO_MATERNO{i},FPTP2003.APELLIDO_MATERNO)&...
        strcmp(PRI_PR.NOMBRE{i},FPTP2003.NOMBRE))+...
        sum(strcmp(PRI_PR.APELLIDO_PATERNO{i},PR2003.APELLIDO_PATERNO)&...
        strcmp(PRI_PR.APELLIDO_MATERNO{i},PR2003.APELLIDO_MATERNO)&...
        strcmp(PRI_PR.NOMBRE{i},PR2003.NOMBRE));
end

% PAN
for i = 1:size(PAN_FPTP,1)
    PAN_FPTP.ran2009(i) = PAN_FPTP.ran2009(i)+...
        sum(strcmp(PAN_FPTP.APELLIDO_PATERNO{i},FPTP2009.APELLIDO_PATERNO)&...
        strcmp(PAN_FPTP.APELLIDO_MATERNO{i},FPTP2009.APELLIDO_MATERNO)&...
        strcmp(PAN_FPTP.NOMBRE{i},FPTP2009.NOMBRE))+...
        sum(strcmp(PAN_FPTP.APELLIDO_PATERNO{i},PR2009.APELLIDO_PATERNO)&...
        strcmp(PAN_FPTP.APELLIDO_MATERNO{i},PR2009.APELLIDO_MATERNO)&...
        strcmp(PAN_FPTP.NOMBRE{i},PR2009.NOMBRE));
    PAN_FPTP.ran2006(i) = PAN_FPTP.ran2006(i)+...
        sum(strcmp(PAN_FPTP.APELLIDO_PATERNO{i},FPTP2006.APELLIDO_PATERNO)&...
        strcmp(PAN_FPTP.APELLIDO_MATERNO{i},FPTP2006.APELLIDO_MATERNO)&...
        strcmp(PAN_FPTP.NOMBRE{i},FPTP2006.NOMBRE))+...
        sum(strcmp(PAN_FPTP.APELLIDO_PATERNO{i},PR2006.APELLIDO_PATERNO)&...
        strcmp(PAN_FPTP.APELLIDO_MATERNO{i},PR2006.APELLIDO_MATERNO)&...
        strcmp(PAN_FPTP.NOMBRE{i},PR2006.NOMBRE))+...
        sum(strcmp(PAN_FPTP.APELLIDO_PATERNO{i},FPTP2006S.APELLIDO_PATERNO)&...
        strcmp(PAN_FPTP.APELLIDO_MATERNO{i},FPTP2006S.APELLIDO_MATERNO)&...
        strcmp(PAN_FPTP.NOMBRE{i},FPTP2006S.NOMBRE))+...
        sum(strcmp(PAN_FPTP.APELLIDO_PATERNO{i},PR2006S.APELLIDO_PATERNO)&...
        strcmp(PAN_FPTP.APELLIDO_MATERNO{i},PR2006S.APELLIDO_MATERNO)&...
        strcmp(PAN_FPTP.NOMBRE{i},PR2006S.NOMBRE));
    PAN_FPTP.ran2003(i) = PAN_FPTP.ran2003(i)+...
        sum(strcmp(PAN_FPTP.APELLIDO_PATERNO{i},FPTP2003.APELLIDO_PATERNO)&...
        strcmp(PAN_FPTP.APELLIDO_MATERNO{i},FPTP2003.APELLIDO_MATERNO)&...
        strcmp(PAN_FPTP.NOMBRE{i},FPTP2003.NOMBRE))+...
        sum(strcmp(PAN_FPTP.APELLIDO_PATERNO{i},PR2003.APELLIDO_PATERNO)&...
        strcmp(PAN_FPTP.APELLIDO_MATERNO{i},PR2003.APELLIDO_MATERNO)&...
        strcmp(PAN_FPTP.NOMBRE{i},PR2003.NOMBRE));
end
for i = 1:size(PAN_PR,1)
    PAN_PR.ran2009(i) = PAN_PR.ran2009(i)+...
        sum(strcmp(PAN_PR.APELLIDO_PATERNO{i},FPTP2009.APELLIDO_PATERNO)&...
        strcmp(PAN_PR.APELLIDO_MATERNO{i},FPTP2009.APELLIDO_MATERNO)&...
        strcmp(PAN_PR.NOMBRE{i},FPTP2009.NOMBRE))+...
        sum(strcmp(PAN_PR.APELLIDO_PATERNO{i},PR2009.APELLIDO_PATERNO)&...
        strcmp(PAN_PR.APELLIDO_MATERNO{i},PR2009.APELLIDO_MATERNO)&...
        strcmp(PAN_PR.NOMBRE{i},PR2009.NOMBRE));
    PAN_PR.ran2006(i) = PAN_PR.ran2006(i)+...
        sum(strcmp(PAN_PR.APELLIDO_PATERNO{i},FPTP2006.APELLIDO_PATERNO)&...
        strcmp(PAN_PR.APELLIDO_MATERNO{i},FPTP2006.APELLIDO_MATERNO)&...
        strcmp(PAN_PR.NOMBRE{i},FPTP2006.NOMBRE))+...
        sum(strcmp(PAN_PR.APELLIDO_PATERNO{i},PR2006.APELLIDO_PATERNO)&...
        strcmp(PAN_PR.APELLIDO_MATERNO{i},PR2006.APELLIDO_MATERNO)&...
        strcmp(PAN_PR.NOMBRE{i},PR2006.NOMBRE))+...
        sum(strcmp(PAN_PR.APELLIDO_PATERNO{i},FPTP2006S.APELLIDO_PATERNO)&...
        strcmp(PAN_PR.APELLIDO_MATERNO{i},FPTP2006S.APELLIDO_MATERNO)&...
        strcmp(PAN_PR.NOMBRE{i},FPTP2006S.NOMBRE))+...
        sum(strcmp(PAN_PR.APELLIDO_PATERNO{i},PR2006S.APELLIDO_PATERNO)&...
        strcmp(PAN_PR.APELLIDO_MATERNO{i},PR2006S.APELLIDO_MATERNO)&...
        strcmp(PAN_PR.NOMBRE{i},PR2006S.NOMBRE));
    PAN_PR.ran2003(i) = PAN_PR.ran2003(i)+...
        sum(strcmp(PAN_PR.APELLIDO_PATERNO{i},FPTP2003.APELLIDO_PATERNO)&...
        strcmp(PAN_PR.APELLIDO_MATERNO{i},FPTP2003.APELLIDO_MATERNO)&...
        strcmp(PAN_PR.NOMBRE{i},FPTP2003.NOMBRE))+...
        sum(strcmp(PAN_PR.APELLIDO_PATERNO{i},PR2003.APELLIDO_PATERNO)&...
        strcmp(PAN_PR.APELLIDO_MATERNO{i},PR2003.APELLIDO_MATERNO)&...
        strcmp(PAN_PR.NOMBRE{i},PR2003.NOMBRE));
end
clear i

tableA3 = table(round(100 * [mean(MP_FPTP.ran2009==1);mean(MP_PR.ran2009==1);...
    mean(NA_FPTP.ran2009==1);mean(NA_PR.ran2009==1);...
    mean(PVEM_FPTP.ran2009(1:N0)==1);mean(PVEM_FPTP.ran2009(N0+1:end)==1);mean(PVEM_PR.ran2009==1);...
    mean(PRI_FPTP.ran2009(1:N0)==1);mean(PRI_FPTP.ran2009(N0+1:end)==1);mean(PRI_PR.ran2009==1);...
    mean(PAN_FPTP.ran2009==1);mean(PAN_PR.ran2009==1)],1),...
    round(100 * [mean(MP_FPTP.ran2006==1);mean(MP_PR.ran2006==1);...
    mean(NA_FPTP.ran2006==1);mean(NA_PR.ran2006==1);...
    mean(PVEM_FPTP.ran2006(1:N0)==1);mean(PVEM_FPTP.ran2006(N0+1:end)==1);mean(PVEM_PR.ran2006==1);...
    mean(PRI_FPTP.ran2006(1:N0)==1);mean(PRI_FPTP.ran2006(N0+1:end)==1);mean(PRI_PR.ran2006==1);...
    mean(PAN_FPTP.ran2006==1);mean(PAN_PR.ran2006==1)],1),...
    round(100 * [mean(MP_FPTP.ran2003==1);mean(MP_PR.ran2003==1);...
    mean(NA_FPTP.ran2003==1);mean(NA_PR.ran2003==1);...
    mean(PVEM_FPTP.ran2003(1:N0)==1);mean(PVEM_FPTP.ran2003(N0+1:end)==1);mean(PVEM_PR.ran2003==1);...
    mean(PRI_FPTP.ran2003(1:N0)==1);mean(PRI_FPTP.ran2003(N0+1:end)==1);mean(PRI_PR.ran2003==1);...
    mean(PAN_FPTP.ran2003==1);mean(PAN_PR.ran2003==1)],1),...
    round(100 * [mean(MP_FPTP.ran2009==1|MP_FPTP.ran2006==1|MP_FPTP.ran2003==1);
    mean(MP_PR.ran2009==1|MP_PR.ran2006==1|MP_PR.ran2003==1);...
    mean(NA_FPTP.ran2009==1|NA_FPTP.ran2006==1|NA_FPTP.ran2003==1);...
    mean(NA_PR.ran2009==1|NA_PR.ran2006==1|NA_PR.ran2003==1);...
    mean(PVEM_FPTP.ran2009(1:N0)==1|PVEM_FPTP.ran2006(1:N0)==1|PVEM_FPTP.ran2003(1:N0)==1);...
    mean(PVEM_FPTP.ran2009(N0+1:end)==1|PVEM_FPTP.ran2006(N0+1:end)==1|PVEM_FPTP.ran2003(N0+1:end)==1);...
    mean(PVEM_PR.ran2009==1|PVEM_PR.ran2006==1|PVEM_PR.ran2003==1);...
    mean(PRI_FPTP.ran2009(1:N0)==1|PRI_FPTP.ran2006(1:N0)==1|PRI_FPTP.ran2003(1:N0)==1);...
    mean(PRI_FPTP.ran2009(N0+1:end)==1|PRI_FPTP.ran2006(N0+1:end)==1|PRI_FPTP.ran2003(N0+1:end)==1);...
    mean(PRI_PR.ran2009==1|PRI_PR.ran2006==1|PRI_PR.ran2003==1);...
    mean(PAN_FPTP.ran2009==1|PAN_FPTP.ran2006==1|PAN_FPTP.ran2003==1);...
    mean(PAN_PR.ran2009==1|PAN_PR.ran2006==1|PAN_PR.ran2003==1)],1),...
    'VariableNames',{'ran2009','ran2006','ran2003','ran2003-2009'},...
    'RowNames',{'MP_FPTP','MP_PR','NA_FPTP','NA_PR','PVEM_FPTPalone','PVEM_FPTPcoalition','PVEM_PR','PRI_FPTPalone','PRI_FPTPcoalition','PRI_PR','PAN_FPTP','PAN_PR',});

disp(' ')
disp('for Table A3:')
disp(tableA3)



